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ABSTRACT 



We present a 3D MHD simulation that qualitatively models the coronal mag- 
netic field evolution associated with the eruptive fiare that occurred on December 
13, 2006 in the emerging 5-sunspot region NOAA 10930 observed by the Hinode 
satellite. The simulation is set up where we drive the emergence of an east-west 
oriented magnetic fiux rope at the lower boundary into a pre-existing coronal 
field constructed from the SOHO/MDI full-disk magnetogram at 20:51:01 UT 
on December 12, 2006. The resulting coronal fiux rope embedded in the ambi- 
ent coronal magnetic field first settles into a stage of quasi-static rise, and then 
undergoes a dynamic eruption, with the leading edge of the fiux rope cavity accel- 
erating to a steady speed of about 830 km/s. The pre-eruption coronal magnetic 
field shows morphology that is in qualitative agreement with that seen in the 
Hinode soft X-ray observation in both the magnetic connectivity as well as the 
development of an inverse-S shaped X-ray sigmoid. We examine the properties 
of the erupting fiux rope and the morphology of the post-reconnection loops, and 
compare them with the observations. 



Introduction 



Coronal mass ejections (CMEs) are large-scale, spontaneous ejections of plasma and 
magnetic fiux from the lower solar corona into interplanetary space and are major drivers 



of space weather near earth (e.g. Hundhausen 1993 Lindsay et al. 1999 Webb et al. 2000). 



CMEs and eruptive fiares are believed to result from a sudden, explosive release of the free 
magnetic energy stored in the previously quasi-equilibrium, twisted/sheared coronal mag- 
netic field (see e.g. reviews by Forbes et al.||2006 Chen|[20lT ). Using idealized constructions, 
both analytical studies and numerical simulations have been carried out to understand the 
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basic underlying magnetic field structures of the eruption precursors, and the physical mech- 



anisms of their sudden eruption (e.g. Mikic and Linker 1994; Antiochos et al. 1999 Forbes 


and Priest||1995 


Lin et al.||1998| |Amari et al.||2000 


Sturrock et al.||2001; Roussev et al.||2003 


Torok and Kliem 2005 2007 Fan and Gibson 2007 Isenberg and Forbes 2007 Fan 2010 


Aulanier et al. 


2010; Demoulin and Aulanier 20 IC 


). Magneto- hydrodynamic (MHD) mod- 



els of observed CME events have also been constructed to determine the actual magnetic 
field evolution and causes for the eruption and the properties of the magnetic ejecta, which 
are critical for determining the geo-effectiveness of the resulting interplanetary coronal mass 



Mikic et al. 


2008 


Titov et al. 


2008 



The eruptive event in active region 10930 on December 13, 2006 produced an X3.4 
flare and a fast, earth-directed CME with an estimated speed of at least 1774 km/s. The 
ICME reached the Earth on 14-15 December 2006, with a strong and prolonged south- 
ward directed magnetic field in the magnetic cloud, causing a major geomagnetic storm 



(e.g. Liu et al. 2008 Kataoka et al. 2009). This event is particularly well observed by 
Hinode for both the coronal evolution as well as the photospheric magnetic field evolution 
over a period of several days preceding, during, and after the eruption. The photospheric 
magnetic field evolution of AR 10930 was characterized by an emerging 5-sunspot with 
a growing positive polarity, which displayed substantial (counter-clockwise) rotation and 
eastward motion as it grew (see e.g. the movie provided at the NOAJ website http://solar^ 



b. nao.ac.jp/news/070321Flare/me_20061208_15arrow_6fps.mpg and see also Min and Chae 



(2009)). This is indicative of the emergence of a twisted magnetic flux rope with the positive 



rotating spot being one of its photospheric footpoints. The total rotation of the positive. 



growing sunspot prior to the onset of the flare is measured to be 240° by Zhang, Li, and 



Song (2007) and 540° by Min and Chae (2009), which gives an estimate of the minimum 



amount of twist that has been transported into the corona in the emerged flux rope. 

Several studies based on non-linear force-free field extrapolations from the photospheric 
vector magnetic field measurement for AR 10930 have been carried out to study the coro- 
nal magnetic field and the associated free magnetic energy before and after the flare (e.g. 



Schrijver et al. 2008 Inoue et al. 2008). In this paper, we present an MHD simulation that 
model the coronal magnetic field evolution associated with the onset of the eruptive flare 
in AR 10930 on December 13 2006. The simulation assumes the emergence of an east-west 
oriented magnetic flux rope into a pre-existing coronal magnetic field constructed based on 
the SOHO/MDI full-disk magnetogram of the photospheric magnetic field at 20:51:01 UT 
on December 12. Our simulated coronal magnetic field first achieves a quasi-equilibrium 
phase during which the coronal flux rope rises quasi-statically as more twisted flux is being 
transported into the corona through a slow flux emergence. The evolution is then followed by 
a dynamic eruption, where the erupting flux rope accelerate to a final steady speed of about 
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830 km/s. The erupting flux rope is found to undergo substantial writhing or rotational 
motion, and the erupting trajectory is non-radial, being deflected southward and eastward 
from the local radial direction of the source region. The coronal magnetic fleld structure 
just prior to the onset of the eruption reproduces qualitatively the observed morphology and 
connectivity of the coronal magnetic fleld, including the formation of an inverse-S shaped 
pre-eruption sigmoid, as seen in the Hinode XRT images. After the onset of the eruption, the 
evolution of the post-reconnection loops and their foot-points resulting from the simulated 
magnetic fleld is also in qualitative agreement with the morphology of the observed X-ray 
post-flare brightening and the evolution of the chromosphere flare ribbons. 

We organize the remainder of the paper as follows. In Section [2| we describe the 
MHD numerical model and how the simulation is set up. In Section [3] we describe the 
resulting evolution of the simulated coronal magnetic fleld and compare with observations. 
We summarize the conclusions and discuss future directions for improving the model in 
Section HI 



2. Model Description 



For the simulation carried out in this study, we solve the following magneto-hydrodynamic 
equations in a spherical domain: 
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In the above v, B, p, p, T, e, e i?, p, 7, G, and Mq denote respectively the velocity field, the 
magnetic field, density, pressure, temperature, the internal energy density, the total energy 
density (internal+kinetic+magnetic), the gas constant, the mean molecular weight, the ratio 
of specific heats, the gravitational constant, and the mass of the Sun. We have assumed 
an ideal polytropic gas with 7 = 1.1 for the corona plasma. The above MHD equations 
are solved numerically without any explicit viscosity, magnetic diffusion, and non-adiabatic 
effects. However numerical dissipations are present, and since we are solving the total energy 
equation in conservative form, the numerical dissipation of kinetic, and magnetic energy is 
effectively being put back into the internal energy. 

The basic numerical schemes we use to solve the above MHD equations are as follows. 
The equations are discretized in spherical domain with r, 6, coordinates using a staggered 



finite-difference scheme (Stone and Norman 1992a), and advanced in time with an explicit. 



second order accurate, two-step predictor-corrector time stepping. A modified, second order 



accurate Lax-Friedrichs scheme similar to that described in Rempel, Schiissler, and Knolker 



(2009, see eq. (A3) in that paper) is applied for evaluating the fiuxes in the continuity 
and energy equations. Compared to the standard second order Lax-Friedrichs scheme, this 
scheme significantly reduces numerical diffusivity for regions of smooth variation, while re- 
taining the same robustness in regions of shocks. The standard second order Lax-Friedrichs 
scheme is used for evaluating the fluxes in the momentum equation. A method of character- 



istics that is upwind in the Alfven waves (Stone and Norman 1992b) is used for evaluating 



the V X B term in the induction equation, and the constrained transport scheme is used to 
ensure V • B = to the machine precision. 

The simulation is set up where we drive the emergence of a part of a twisted magnetic 
torus at the lower boundary into a pre-existing coronal potential field, constructed based 
on the MDI full-disk magnetogram from 20:51:01 UT on December 12, 2006 (Figure [l^). 
First, from the full-disk MDI magnetogram, a region centered on the 5-spot (the white box 
in Figure [1^), with an latitudinal extent of 30° and a longitudinal extent of 45° is extracted 
as the lower boundary of the spherical simulation domain. In terms of the the simulation 
coordinates, the domain spans r G [Rq,Q.25Rq], 6 e [75°, 105°], (f) G [-22.5°, 22.5°], with 
the center of its lower boundary: 9 = 90° and = 0°, corresponding to the center of 
the white-boxed area in Figure [l^. This domain is resolved by a grid of 512 x 352 x 528, 
with 512 grid points in r, 352 grid points in 6, and 528 grid points in 0. The grid is 
uniform in the 6 and directions but non-uniform in r, with a uniform grid spacing of 
dr = 1.028 Mm in the range of r = Rq to about 1.6i?0 and a geometrically increasing grid 
spacing above 1.6i?0, reaching about dr = 173.4 Mm at the outer boundary. We assume 
perfectly conducting walls for the side boundaries, and for the outer boundary we use a 
simple outward extrapolating boundary condition that allows plasma and magnetic field to 
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flow through. The lower boundary region extracted from the MDI full disk magnetogram (as 
viewed straight-on) is shown in Figure where we simply take the interpolated line-of-sight 
flux density from the full-disk magnetogram and assume that the magnetic field is normal 
to the surface to obtain the Br shown in the Figure. The region contains roughly all the flux 
of the (5-spot and the surrounding pores and plages, to which some of the flux of the (5-spot 
is connected. The peak field strength in the region is about 3000 G. A smoothing using a 
Gaussian filter is carried out on the lower boundary region until the peak field strength is 
reduced to about 200 G. This is necessary since the simulation domain corresponds to the 
corona, with the lower boundary density assumed to be that of the base of the corona, and 
thus a significant reduction of the field strength from that measured on the photosphere is 
needed to avoid unreasonably high Alfven speeds, which would put too severe a limit on 
the time step of numerical integration. After the smoothing, the magnetic flux in a central 
area, which roughly encompasses the region of the observed flux emergence (including the 
rotating, positive sunspot) is zeroed out (see Figure [!}:;) to be the area where the emergence 
of an idealized, twisted magnetic torus is driven on the lower boundary. The potential field 
constructed from this lower boundary normal flux distribution in Figure ^ is assumed to be 
the initial coronal magnetic field for our simulation, which is shown in Figure |2j We zero out 
the normal flux in the area for driving the flux emergence so that we can specify analytically 
the subsurface emergence structure in a field free region without the complication of the 
subsurface extension of a pre-existing flux in the same area. 

The initial atmosphere in the domain is assumed to be a static polytropic gas: 



where po = 8.365 x 10^^^ g cm~^, and po = 0.152 dyne cm~^ are respectively the density 
and pressure at the lower boundary of the coronal domain, and the corresponding assumed 
temperature at the lower boundary is 1.1 MK. The initial magnetic field in the domain 
is potential, and thus does not exert any forcing on the atmosphere which is in hydrostatic 



a vertical line rooted in the peak Bj. of the main pre-existing negative polarity spot. For 
the initial state constructed, the peak Alfven speed is about 24 Mm/s, and the sound speed 
is 141 km/s at the bottom and gradually declines with height. In most of the simulation 
domain, the Alfven speed is significantly greater than the sound speed. 

At the lower boundary (at r = -R0), we impose (kinematically) the emergence of a 
twisted torus Btube by specifying a time dependent transverse electric field E±\r=RQ that 




(9) 




(10) 



equilibrium. Figure |3] shows the height profiles of the Alfven speed and the sound speed along 
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corresponds to the upward advection of the torus with a velocity Vf 



-^Vhsc X Btnhc{RQ,0, (f),t) ] X r 



'11^ 



The magnetic field Btube used for specifying E±\r=RQ is an axisymmetric torus defined in 
its own local spherical polar coordinate system (r', 6', (f)') whose polar axis is the symmetric 
axis of the torus. In the sun-centered simulation spherical coordinate system, the origin of 
the (r', 9', </)') system is located at r = = {rc,6c,(f)c), and its polar axis (the symmetric 
axis of the torus) is in the plane of the 9 and vectors at position Vc and tilted from the —9 
direction clockwise (towards the direction) by an angle 6. In the (r', 9', 0') system, 

Btube = V X ( 4^^0' 1 + Bs,(r', 9'U\ (12) 
r' sm 9' 



where 



A(r',e') = \,a^B.(l-^^] , (13) 



BAr'J')^^(l-^). (14) 

r sm 6'' \ J 

In the above, a is the minor radius of the torus, uj = (r'^ + R'^ — 2r'i?' sin ^')^/^ is the 
distance to the curved axis of the torus, where R' is the major radius of the torus, q denotes 
the angular amount (in rad) of field line rotation about the axis over a distance a along 
the axis, and Bta/R' gives the field strength at the curved axis of the torus. The magnetic 
field Btube is truncated to zero outside of the fiux surface whose distance to the torus axis 
isvu = a. We use a = O.OSSi?©, R' = OMSRq, q/a = -0.0308 rad Mm-\ Bta/R' = 111 G. 
The torus center is assumed to be initially located at Vc = (rc = 0.902/?©, 9c = 90°, 0c = 0°), 
and the tilt of the torus 6 = 0. Thus the torus is initially entirely below the lower boundary 
and is in the azimuthal plane. For specifying E±\r=RQ, we assume that the torus moves 
bodily towards the lower boundary at a velocity Vfisc = "^^riserc, where v^se is described 
later. The imposed velocity field at the lower boundary is a constant Vfise in the area where 
the emerging torus intersects the lower boundary and zero in the rest of the area. The 
resulting normal flux distribution on the lower boundary after the imposed emergence has 
stopped is shown in Figure [T]i. In it an east-west oriented bipolar pair has emerged, where 
the positive spot represents the emerging, rotating positive sunspot at the south edge of 
the dominant negative spot in Figure [T]d, and the negative spot corresponds to the flux in 
the fragmented pores and plages to the west of the rotating positive sunspot in Figure [TJd. 



Observational study by Min and Chae (2009) found that the minor, fragmented pores of 



negative polarity emerged and moved westward while the positive rotating sunspot moved 
eastward, suggesting that they are the counterpart to which the positive rotating sunspot 
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is at least partly connected to (see Figure 2 in Min and Chae (2009)). This is one of the 
reasons that we model the coronal magnetic field in this study with the emergence of an 
east-west oriented twisted fiux rope. After the emergence is stopped, the transverse electric 
field on the lower boundary (eq. [iT]) is set to zero and the magnetic field is line-tied at the 
lower boundary. At the end of the emergence, the peak normal field strength in the emerged 
bipolar region on the lower boundary reaches 121 G, compared to the 178 G peak normal 
field strength in the dominant negative pre-existing spot in the initial lower boundary field. 
Due to the substantial smoothing of the observed normal magnetic fiux density, the total 
unsigned fiux on the lower boundary of our simulation is only about 30% of that on the 
photosphere in the boxed area shown in Figure [T^. However, the ratio of the emerged fiux 
(in the fiux rope) over the total fiux on the lower boundary, ~ 10%, for our simulation is 
about the same as the ratio of the observed emerged fiux (in the positive rotating sunspot) 
over the total fiux in the boxed area in Figure [l^. 

Note that although the coronal temperature and density are used at the lower boundary, 
the dynamic property of the lower boundary refiects the property of the photosphere. The 
lower boundary is assumed to be "infinitely heavy" such that the magnetic stress exerted on 
it from the corona does not result in any motion of the field line foot-points (field anchoring 
or line-tying) and that the lower boundary evolves in a prescribed way by a kinematically 
imposed fiux emergence associated with the upward advection of a twisted fiux rope. Thus 
dynamically the lower boundary is meant to approximate the photosphere, which can sup- 
port cross-field currents and the resulting magnetic stresses. However the thermodynamic 
conditions of the corona (instead of the photosphere) are used for the lower boundary so that 
(1) we do not have to resolve the small (about 150 km) photospheric pressure scale height 
in a simulation of the large scale coronal evolution of a CME (size scale on the order of a 
solar radius), and (2) we avoid solving the complex energy transport associated with coronal 
heating, radiative cooling, and thermal conduction, which would be required if we were to 
include the thermodynamics of the photosphere-chromosphere-corona system in the simula- 
tion. Here for modeling the large scale, magnetically dominated dynamic evolution of the 
CME initiation, we greatly simplify the thermodynamics (assuming an ideal polytropic gas 
for the coronal plasma throughout the domain), and focus on the magnetic field evolution of 
the corona in response to the imposed fiux emergence and field-line anchoring representative 
of the heavy photospheric lower boundary. 

In the remainder of the paper, quantities are expressed in the following units unless 
otherwise specified: Rq = 6.96 x 10^° cm, po = 8.365 x 10"^^g/cm^, Bq = 20 G, Vao = 
Bol ^/'i'Kpo = 1951 km/s, r^o = -Ro/fao = 356.8 s, as units for length, density, magnetic 
field, velocity and time respectively. Due to the large peak Alfven speed (~ l2vao ~ 24, 000 
km/s) in the domain (see Figure [3]), we initially drive the emergence of the twisted torus 
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through the lower boundary at a fairly high speed over a period oft = 0tot = 1.2 with 
Wrise = O.OSfao ~ 98km/s, which is just under the sound speed at the lower boundary but 
significantly slower than the Alfven speed. In this way we build up the pre-eruption coronal 
magnetic field approximately quasi-statically and yet fast enough to minimize numerical 
diffusion. After t = 1.2, we significantly reduce the driving speed of the flux emergence at 
the lower boundary to Wrise = O.Olfao and thus allow the coronal magnetic field to evolve 
quasi-statically until it erupts dynamically. 



3. Results 

Figures |4] and [s] show snapshots of the 3D coronal magnetic field evolution (as viewed 
from 2 different perspectives) after the initial stage of relatively fast emergence has ended 
at if: = 1.2, and the speed for driving the flux emergence at the lower boundary has been 
reduced to v^se = O.OlVao- The view shown in Figure |4] corresponds to the observation 
perspective at the time of the flare, for which the center of the emerging region (also the 
center of the simulation lower boundary) is located at 7.1°S and 24°W from the solar disk 
center (or the line-of-sight). GIF movies for the evolution shown in Figure |4] and Figure |5] 
are available in the electronic version of the paper. We see that the emerged coronal flux 
rope settles into a quasi-static rise phase and then undergoes a dynamic eruption. Figure |6] 
shows the evolution of the rise velocity Vr measured at the apex of the tracked axis of the 
emerged flux rope (triangle points), and also measured at the leading edge of the flux rope 
(crosses). After the emergence is slowed down at t = 1.2, the rise velocity at the apex of 
the flux rope axis slows down, and undergoes some small oscillations as the flux rope settles 
into a quasi-static rise. The quasi-static rise phase extends from about t = 1.2 until about 
t = 2.5, over a time period of 1.3, long compared to the dynamic time scale of ~ 0.1 for the 
estimated Alfven crossing time of the flux rope. At about t = 2.5, the flux rope axis starts 
to accelerate significantly and a dynamic eruption ensues. The flux emergence is stopped 
at t = 2.8, after which the flux rope continues to accelerate outward. We are able to follow 
the acceleration of the axial field line up to = 0.54 = 1050 km/s at t = 3., when the 
axial field line undergoes a reconnection and we are subsequently unable to track it. Figure 
[6] also shows Vr measured at the leading edge of the low density cavity (as shown in Figure 
[7|, corresponding to the expanding flux rope. We find that by the time of about t = 3.2, 
a shock front followed by a condensed sheath has formed ahead of the flux rope cavity (see 
Figure [T] at t = 3.25), and the Vr measured at the front edge of the cavity (or the inner edge 
of the sheath) reaches a steady speed of about 0.425 or 830 km/s (see crosses in Figure [6]). 

When the flux rope begins significant acceleration (at t ^ 2.5), the decay index n = 
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ciln |Bp|/(iln/i which describes the rate of dechne of the corresponding potential field Bp 
with height h is found to be n ^ 1.2 at the apex of the flux rope axis, and n ^ 1.4 at the 
apex of the flux rope cavity. These values are smaller than the critical value of Ucnt = 1-5 



for the onset of the torus instability for a circular current ring (Bateman 1978 Kliem and 



Torok 2006 Demoulin and Aulanier 2010), although there is a range of variability for the 



critical value ricnt, which can be as low as 1, depending on the shape of the current channel 



of the flux rope (e.g. Demoulin and Aulanier 2010). For a 3D anchored flux rope, as is the 
case here, it is difficult to obtain an analytical determination of ricrit for the instability or loss 



of equilibrium of the flux rope (Isenberg and Forbes 2007). The exact critical point for the 



onset of the torus instability would depend on the detailed 3D magnetic field configuration. 
On the other hand, a substantial amount of twist has been transported into the corona at 
the onset of eruption. At t = 2.5, the self-helicity of the emerged fiux rope reaches about 
— 1.02$^Qpg, where $rope is the total magnetic fiux in the rope, corresponding to field lines in 
the fiux rope winding about the central axis by about 1.02 rotations between the anchored 
foot points. This suggests the possible development of the helical kink instability of the 



fiux rope (e.g. Hood and Priest 1981 Torok and Kliem 2005 Fan and Gibson 2007). The 
erupting fiux rope is found to undergo substantial writhing or kinking motion as can be seen 
in the sequences of images (also the movies in the electronic version) in Figures 111 and IS] 



We also find that the trajectory for the eruption of the fiux rope is not radial because 
of the ambient coronal magnetic field: the erupting fiux rope is defiected southward and 
eastward from the local radial direction (see Figures |4| [sj and [7| and the associated movies). 
Using the apex location of the erupting fiux rope cavity at t = 3.25 (Figure [T]), we find 
that the erupting trajectory at that time is defiected by 2.3° southward and 1.3° eastward 
from the local radial direction at the center of fiux emergence, and further defiection of the 
trajectory continues with time. Since the local radial direction at the center of the fiux 
emergence corresponds to 7.1°S and 24° W from the solar disk center (or the line-of-sight), 
the defiection during the eruption is sending the fiux rope towards the line-of-sight in the 
east-west direction, but further southward away from the line-of-sight in the north-south 
direction. This is consistent with the observed halo of the CME seen in LASCO C2 and 



C3 coronagraphs (Figure 2 in Kataoka et al. (2009) and Figure 1 in Ravindra and Howard 



(2010)), where the north-south and east-west asymmetries of the halo distribution indicate 
that the direction of ejection is more southward and less westward than what would have 
been expected for a radial ejection from the location of the source region on the solar disk. 

Figure [s] shows the coronal magnetic field as viewed from the side (panels a and b) and 
viewed from the observing perspective (panels c and d) just before the onset of eruption at 
t = 2.45, compared with the Hinode XRT image of the region (panel e) just before the fiare. 
We see that the morphology of the coronal magnetic field and its connectivity are very similar 
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to those shown in the X-ray image. To understand the nature of the bright X-ray sigmoid 
in the image, we have identified the region of significant magnetic energy dissipation and 
heating in the simulated magnetic field using both the electric current density J = |V x B| 
and the increase of entropy AS = C^Aln(p/p^). As pointed out in Section [2| since we are 
solving the total energy equation in conservative form, numerical dissipation of magnetic 
energy and kinetic energy due to the formation of current sheets and other sharp gradients 
is being implicitly put back into the thermal energy of the plasma, resulting in an increase 
of the entropy. We have identified regions where there is significant entropy increase with 
AS/Ct, > 1.15 and also high electric current density concentration with J/B > l/l where 
/ = 10 times the grid size. Such regions are outlined by the orange iso-surfaces in panels (a) 
and (c) of Figure |8| and they appear as an inverse-S shaped layer (as viewed from the top), 
which likely corresponds to the formation of an electric current sheet underlying the anchored 
flux rope (e.g Titov and Demoulin|1999 Low and Berger [2003 [Gibson et al.||2006 ). We have 



also plotted field lines (purple field lines shown in panels b and d) going through the region 
of the current layer, which are preferentially heated and are expected to brighten throughout 
their lengths (due to the high heat conduction along the field lines) in soft-X ray, producing 
the central dominant X-ray sigmoid seen in the Hinode XST image (panel e). Thus our quasi- 
equilibrium coronal magnetic field resulting from the emergence of a nearly east-west oriented 
magnetic fiux rope could reproduce the observed overall morphology and connectivity of the 
coronal magnetic field, including the presence of the observed pre-eruption X-ray sigmoid. 
We find that both J/B as well as A^* peak along the "left elbow" portion of the current 
layer, where the positive polarity fiux of the emerged fiux rope comes in contact with the 
flux of the dominant pre-existing negative polarity sunspot, consistent with the brightness 
distribution along the observed X-ray sigmoid (panel e of Figure [s]) . Reconnections in this 
part of the current layer cause some of the flux in the emerged flux rope to become connected 
with the major negative sunspot (see the green fleld lines connecting between the dominant 
negative spot and the emerging positive spot in panel (d) of Figure [s]). We have also done 
a few simulations where we varied the tilt of the emerging flux rope, and found that to 
reproduce the observed orientation of the sigmoid, the emerging flux rope needs to be nearly 
east-west oriented. 

With the onset of the eruptive flare, the soft-X ray observation flrst shows a transient 
brightening of the sigmoid, and subsequently the emission is completely dominated by the 
brightness of the post-flare loops (see panels (a)(c)(e) of Figure |9]). In the simulated coro- 
nal magnetic fleld, we flnd that the current density in the inverse-S shaped current layer 
intensifles as the flux rope begins to erupt. We can deduce qualitatively the evolution of 
the post-reconnection (or post-flare) loops from our modeled magnetic fleld evolution. We 



traced fleld lines (see the red fleld lines in panels (b) (e) (h) of Figure 10 and panels (b) (d) (f ) 
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of Figure 11) whose apexes are located in the layer of the most intense current density and 
heating. These field lines are the ones who have just reconnected at their apexes and would 
slingshot downwards, corresponding to the downward collapsing post-flare loops. The layer 
of the most intense current density and heating, as outlined by the orange iso-surfaces in 
panels (a)(d)(g) of Figure 10 and panels (a)(c)(e), is identified as where J / B > l/l with 
/ = 5 times the grid size, and where AS/Cy > 2.3. This most intense current layer is 
found to rise upward with the eruption of the flux rope. The associated post-reconnection 
field lines are initially low lying and form a narrow sigmoid shaped bundle as can be seen 
in Figures [lO|b) and[TT|^b). With time, the post-reconnection loops broaden and rise up, 
showing cusped apexes (Figures [lO|^e)(h) and Figures [Il|^d)(f)). The morphology of the 
post-reconnection loops, which transition from an initially narrow low-lying sigmoid bundle 
to a broad, sigmoid-shaped row of loops with cusped apexes is in qualitative agreement with 
the observed evolution of the post-flare X-ray brightening shown in Figures |9](a)(c)(e). 



The foot points of the post-reconnection loops (panels (c)(f)(i) of Figures 10) can be 
compared qualitatively with the evolution of the flare ribbons in the lower solar atmosphere 
as shown in the Hinode SOT observation (panels (b)(d)(f) of Figure |9]). The ribbon cor- 
responding to the positive polarity foot points (orange ribbon in Figures [IO|(c)(f)(i)) of 
the post-reconnection loops is found to sweep southward across the newly emerged posi- 
tive polarity spot, similar to the apparent movement of the positive polarity ribbon seen 
the observation (panels (b) (d) (f ) of Figure |9]) in relation to the observed positive emerging 
spot. For the ribbon corresponding to the negative polarity footvpoints (the yellow ribbon 
in Figures [lO|^c)(f)(i)), its eastern portion is found to extend and sweep northward into the 
dominant pre-existing negative spot, while its western hook-shaped portion is found to sweep 
northward across the newly emerged negative spot. Similarly in the SOT observation (panels 
(b)(d)(f) of Fig [9]), for the negative polarity ribbon, the eastern portion sweeps northward 
into the dominant negative sunspot, while its western, upward curved hook-shaped portion is 
found to sweep northward across the minor, fragmented negative pores which have emerged 
to the west of the main 5-sunspot. The modeled ribbons based on the footpoints certainly 
differ in many ways in their shape and extent compared to the observed flare ribbons. But 
they capture some key qualitative features in the observed motions of the flare ribbons in 
relation to the photospheric magnetic flux concentrations. 



4. Discussions 



We have presented an MHD model that qualitatively describes the coronal magnetic field 
evolution of the eruptive flare in AR 10930 on December 13, 2006. The model assumes the 
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emergence of an east-west oriented magnetic flux rope into a pre-existing coronal magnetic 
field constructed based on the MDI full-disk magnetogram of the photospheric magnetic field 
at 20:51:01 UT on December 12. As described in Section [2| a substantial smoothing of the 
observed photospheric magnetic flux density from the MDI magnetogram is carried out such 
that the peak field strength on the lower boundary is reduced from ~ 3000 G to ~ 200 G 
to avoid the extremely high Alfven speed that would put too severe a limit on the time step 
of numerical integration. The imposed flux emergence at the lower boundary of an idealized 
subsurface magnetic torus produces a flux emergence pattern on the lower boundary that 
is only qualitatively representative of the observed flux emergence pattern (compare Figure 
[T|d and Figure [l|l). In the model, the emerging bipolar pair on the lower boundary is more 
symmetric, more spread-out in spatial extent, and both polarities are transporting left- 
handed twist (or injecting negative helicity flux) into the corona at the same rate. Whereas 
in the observation, the positive emerging sunspot is coherent and clearly shows a counter- 
clockwise twisting motion, indicating an injection of negative helicity flux into the corona. 



while its counterpart to the west is in the form of fragmented pores (e.g. Min and Chae 



2009). However a quantitative measurement by Park et al. (2010) using MDI magnetograms 
also found a significant negative helicity flux associated with these fragmented pores as well 
(see Figure 4 in their paper). In the simulation, the self-helicity of the emerged portion 
of the flux rope in the corona at the end of the imposed flux emergence (at t = 2.8) is 
-f^ropc ~ 1.07$^, where $ is the normal flux in each polarity of the emerged bipolar region 
on the lower boundary. This is a measure of the internal twist in the emerged flux rope and 
it corresponds to field lines twisting about the axis by about 1.07 winds (or 385° rotation) 
between the two anchored ends in the emerged flux rope. On the other hand, the total 
relative magnetic helicity i^tot that has been transported into the corona by the imposed 
flux emergence is found to be i^tot ~ 3.02$^, which is the sum of both the self-helicity of 
the emerged flux rope H^^pc as well as the mutual helicity between the emerged flux and 
the pre-existing coronal magnetic field. The observed amount of rotation of the positive 



emerging sunspot, ranging from 240° (Zhang, Li, and Song 2007) to 540° (Min and Chae 
2009), gives an estimate of (ifrope/'^'^) x 360° for the emerged flux rope, which is about 385° 



in the simulation and is thus within the range of the observed values. 

After an initial phase where we drive the emergence of the twisted torus at a fairly large 
(but still significantly sub-Alfvenic) speed to quickly build up the pre-eruption field, we slow 
down the emergence and the coronal magnetic field settles into a quasi-equilibrium phase, 
during which the coronal flux rope rises quasi-statically as more twist is being transported 
slowly into the corona through continued flux emergence. This phase is followed by a dynamic 
eruption phase where the coronal flux rope accelerates in the dynamic time scale to a steady 
speed of about 830 km/s. Due to the substantial twist (greater than 1 full wind of field line 
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twist) that has been transported into the corona at the onset of the eruption, the erupting 
flux rope is found to undergo substantial writhing motion. The erupting flux rope underwent 
a counter-clockwise rotation that exceeded 90° by the time the front of the flux rope cavity 
reached 1.4i?0. We also find that the initial trajectory of the erupting flux rope is not 
radial, but is deflected southward and eastward from the local radial direction due to the 
ambient coronal magnetic field. Since the initial coronal flux rope is located at 7.1°S and 
24°W from the solar disk center, the deflection is sending the erupting flux rope towards the 
line-of-sight in the east-west direction, but further away from the line-of-sight in the north- 
south direction, consistent with the observed halo of the CME seen in LASCO C2 and C3 
coronagraphs, where the halo's north-south (east-west) asymmetry appears larger (smaller) 
than would have been expected from a radial eruption of the flux rope from its location on 
the solar disk. However, due to the relatively restrictive domain width in 9 (30°) and 
(45°) in our current simulation, the side wall boundary in the south begins to significantly 
constrain the further southward deflection and expansion of the flux rope by the time the top 
of the flux rope cavity reaches about 1.4i?0. Thus, we are not able to accurately determine 
the subsequent trajectory change or the continued writhing of the erupting flux rope beyond 
this point. A larger simulation with a significantly greater domain size in 9 and 0, that still 
adequately resolves the coronal magnetic field in the source region, will be carried out in a 
subsequent study to determine the later properties of the flux rope ejecta. 

The restrictive domain size may also play a role in the significantly lower steady speed 
of 830 km/s reached by the erupting flux rope in the simulation, compared to the observed 



value of at least 1780 km/s for the speed of the CME (e.g. Ravindra and Howard 2010). It 
has been shown that the rate of spatial decline of the ambient potential magnetic field with 
height is both a critical condition for the onset of the torus instability of the coronal flux rope 
( Batemaii||1978 Kliem and Torok~||2006 Isenberg and Forbes||2007 Fan||20ld Demoulin and 



Aulanier 2010 e.g.) as well as an important factor in determining the acceleration and the 



final speed of the CMEs ( Torok and Kliem [2007 ). Even for a sufficiently twisted coronal flux 
rope that is unstable to the helical kink instability, the spatial decline rate of the ambient 
potential field is found to determine whether the non-linear evolution of the kink instability 
leads to a confined eruption (with the flux rope settles into a new kinked equilibrium) or 
an ejection of the flux rope (Torok and Kliem 2005). The simulation in this paper has 



assumed perfect conducting walls for the side boundaries where the field lines are parallel to 
the walls. Thus widening the simulation domain would result in a more rapid expansion and 
hence a more steep decline of the ambient potential field with height. This would result in a 
greater acceleration and a faster final speed for the CME based on the results from previous 



investigations by Torok and Kliem (2005, 2007). It may be difficult to distinguish whether 



the torus or the kink instability initially triggers the eruption given the complex 3D coronal 
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magnetic field, but the final speed of the CME would be strongly affected by the spatial 
decline rate of the ambient potential field for either cases. The substantial smoothing of the 
lower boundary magnetic field to reduce the peak Alfven speed is also a major reason for 
the low final speed of the erupting flux rope in the current simulation. 

Nevertheless, our simulated coronal magnetic field evolution is found to reproduce sev- 
eral key features of the eruptive flare observed by Hinode. The pre-eruption coronal field 
during the quasi-static phase reproduces the observed overall morphology and connectivity 
of the coronal magnetic field, including the presence of the pre-eruption X-ray sigmoid seen 
in the Hinode XRT images. The presence of the pre-eruption sigmoid in our model is caused 
by the preferential heating of an inverse-S shaped flux bundle in the flux rope by the forma- 
tion of an inverse-S shaped current sheet underlying the flux rope. Our simulations suggest 
that the emerging flux rope needs to be nearly east-west oriented in order to reproduce the 



observed orientation of the X-ray sigmoid. This is consistent with the suggestion by (Min 



and Chae 2009 ) that the counterpart of the emerging, rotating positive sunspot is the mi- 
nor negative pores to the west of the emerging sunspot (rather than the dominant negative 
sunspot). After the onset of the eruption, the morphology of the post-fiare loops deduced 
from the simulated field show a transition from an initial narrow, low-lying sigmoid bundle 
to a broad, sigmoid-shaped row of loops with cusped apexes, in qualitative agreement with 
the evolution of the post-fiare X-ray brightening observed by XRT of Hinode. The apparent 
motions of the foot points of the post-fiare loops in relation to the lower boundary magnetic 
flux concentrations are also in qualitative agreement with the evolution of the chromospheric 
flare ribbons observed by Hinode SOT. These agreements suggest that our simulated coronal 
magnetic field produced by the emergence of an east-west oriented twisted flux rope, with the 
positive emerging flux "butting against" the southern edge of the dominant pre-existing neg- 
ative sunspot, captures the gross structure of the actual magnetic field evolution associated 
with the eruptive flare. To improve quantitative agreement, a more accurate determination 



of the lower boundary electric field (Fisher et al. 2011) that more closely reproduces the 



observed flux emergence pattern on the lower boundary is needed. 
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Fig. 1. — (a) MDI full-disk magnetogram taken at 20:51:01 UT on December 12, 2006. The 
surface area enclosed by the white box corresponds to the lower boundary surface of the 
simulation domain, (b) Br on the lower boundary region as viewed straight-on from the 
center of the region, (c) Br on the lower boundary after applying a Gaussian smoothing and 
with the field in a central region being zeroed out for imposing the emergence of a twisted 
magnetic flux rope. The potential field extrapolated from the Br shown here is the assumed 
initial field in the simulation domain (see Fig. |2] below) . (d) Br on the lower boundary at 
the end of emergence of the twisted flux rope. 
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Fig. 2. — Selected field lines of the initial potential magnetic field for the simulation as 
viewed from two different perspectives. 
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Fig. 3. — The Alfven speed and the sound speed as a function of raidial distance along 
a vertical line rooted in the peak of the main pre-existing negative spot on the lower 
boundary 



Fig. 4. — Snapshots of the 3D coronal magnetic field evolution after t = 1.2 (when the initial 
fast emergence phase has ended), showing the quasi-static rise and then the eruption of the 
coronal fiux rope. The field is viewed from the observation perspective at the time of the 
observed fiare, where the center of the emerging fiux rope is located at 7.1°S and 24°W from 
the solar disk center (or the line-of-sight). A .gif movie of the above evolution is available in 
the electronic version of the paper 




Fig. 5. — Same as Figure |4] except viewed from a different perspective. A .gif movie of the 
evolution is available in the electronic version of the paper 




Fig. 6. — The evolution of the radial velocity measured at the apex of the tracked axial field 
hue of the emerged flux rope (triangles) and also measured at the leading edge of the flux 
rope cavity (crosses). 
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Fig. 7. — A meridional slice of density through the middle of the erupting flux rope. It shows 
a low density cavity corresponding to the expanding flux rope. A shock front has formed 
with a dense sheath compressed between the flux rope and the shock front. 
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(a) (b) 




Fig. 8. — Coronal magnetic field as viewed from the side (a and b) and viewed from the 
observing perspective (c and d) just before the onset of eruption at t = 2.45, compared with 
the Hinode XRT image of the region (e) just before the flare. The orange surfaces are the 
iso-surfaces where J / B = l/l with / = 10 times the grid size and where AS/C^ > 1.15. They 
outline the region of strong electric current layers. The purple field lines are the field lines 
that go through the points in the current layer. 
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Fig. 9. — Hinode XRT images of the post-flare brightening (left column images), and the 
corresponding Hinode SOT snapshots in Ca II line showing the chromosphere flare ribbons 
(right column images). The orange (yellow) arrows indicate the extent of the flare ribbon in 
the positive (negative) polarity. 



t = 3.0 t = 3.0 t = 3.0 




Fig. 10. — Left column images: iso-surfaces where J/B = l/l with / = 5 times the grid 
size and where AS/C^ > 2.3, outhning the the most intensely heated portion of the current 
layer. Middle column images: sampled field lines whose apexes are in the intense current layer 
outlined by iso-surfaces in the left column images, corresponding to the post-reconnection 
loops. Right column images: foot points of the post-flare loops shown in the middle column 
images. The gray scale images in all panels show the normal magnetic field at the lower 
boundary of the simulation domain 
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Fig. 11. — Same as the left and middle columns of Figure 10 but viewed from a different 
perspective. 



